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The wave functions pertaining to a bilinearly stratified medium are presented in this 
paper. Solutions to the equation V*[«(r)v^] = have been investigated for several profiles 



of *e(r). An analysis is given to the equation v[/c-V^l= — p/e which arises from 
formulation of the quasi-static electric field in a homogeneous anisotropic medium. 
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1. Introduction 

It is known that the core of all electromagnetic 
problems involving* a spherically stratified medium 
consists of two second order linear differential 
equations of the form: 



d 2 S n 
dr 2 " 



.^.S&H^Jflf. 







(1) 



and 



1 ; dr LkW 



dT„ 
dr 



] + [ft(r)-^tl)]r„=0 (2) 



where 



k(/') = 



e(r). 



Based upon the two sets of scalar functions S n and 
T n , four sets of vector wave functions may be 
constructed to describe the electromagnetic field 
within such a medium. The formal solutions to the 
spherical Luneberg lens [Tai, 1958], the cylindrical 
Luneberg lens [Tai, 1956], and the Maxwell fish eyes 
[Tai, 1958], have already been reported based upon 
this method. We shall present here the solutions to 
(1) and (2) corresponding to a bilinear distribution 
of k(t) that contains the stratification of the conical 
lens of Luneberg [1944] as a special case. 

In studying the scattering of electromagnetic 
waves by small spherically stratified particles, it is 
not necessary to find the complete wave solution. 
By means of the Kayleigh method, or its extension 
[Tai, 1952], it is sufficient to investigate the solution 
to the potential function yf/ which satisfies the 
differential equation 
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y[/c(/-)W] = 0. 



(3) 



We shall summarize here those cases in which the 
wave functions have already been found for the same 
profile of /c(r). 

Finally, a brief discussion will be given to the 
solution of the equation 



V '[k-vH= —pfco> 



(4) 



where k denotes a homogeneous dielectric tensor. 

Equation (4) arises from a consideration of the 

quasi-static solution to the electric field in an 
anisotropic medium. 

2. Wave Functions Pertaining to a 
Bilinearly Stratified Medium 

A bilinearly stratified medium is defined as one 
with the dielectric "constant" varying according to 
the following relation: 



By writing 
(5) becomes 



K 



(5) 



•*-*- CO 



k(t)=K„ 



r+ 



(S 



r 2 



r+r 2 



(6) 



where K and K m denote, respectively, the values of 
k(v) defined at r=0 and r= c° . At r=r 2 , the value of 
k(v 2 ) is equal to the mean of K and K^. One may, 
therefore, choose the proper values of K , K a1 and r 2 
to simulate a monotonically increasing or decreasing 
function of k(t) which has a finite asymptotic value 
at r= oo . 
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For convenience, the independent variable r in (1) 
and (2) will be changed to p by letting p—kr, then 
(5) becomes 



K( P )=K a 



p-rp2 



(7) 



where pi=kr, and p 2 = kr 2 . By substituting (7) into 
(1), and letting 



S n =(p + P2 ) P n+1 Un, 



(8) 



one finds that the function U n satisfies the following 
differential equation: 



where 



2n + \ 



P + Pi/ \ P P + P2/ 



(9) 



a=K 03 , a -- 



2n+l 
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ai=K a (pi — p 2 )' 



271+1 

; 

P2 



and the prime on U n in (9) denotes the derivative 
with respect to p. The function U n has two regular 
singularities and one irregular singularity at infinity 
of the second species. Adopting the notation of Ince 
[1944], one would designate this by [0,2, 1 2 ]. The 
same type of function occurs in the theory of the 
spherical Luneberg lens [Tai, 1958]. However, in 
that case it was the T n function which is of this type. 
The S n function in the theory of the Luneberg lens 
is expressible in terms of the confluent hypergeo- 
metric function. 

By substituting (7) into (2) and letting 



l V n , 



(10) 



one finds that V n satisfies the following equation: 
2n+l 1 



K'+( 2 



where 



+ 



P + Pl ' P + P2/ 

+(b+ bo > bl ' 



P + Pl 



-P2/ 



V n =Q, (11) 



\P2 PJ 



b= K a> ,b a =i2, 



6 1= ^±1, .*=*„(*-»). 2n+1 



Pi 



P2 



The function V n has one more regular singular at 
p=—p 2 . It is of the type [0,3, 1 2 ]. 

When one puts r 2 ==0 in (7), the distribution of k(t) 
becomes the same as the profile of the conical lens of 
Luneberg [1944]. The U n and V n functions then re- 
duce, respectively, to the types [0,1, 1 2 ] and [0,2,1 2 ]. 
The type [0,1, 1 2 ] is simply the confluent hypergeo- 



metric function. It is interesting to observe that the 
same two types of functions are involved in the 
electromagnetic theory of the spherical Luneberg 
lens. 

3. Some Potential Functions 

In either purely electrostatic problems dealing with 
a spherically stratified medium or scattering problems 
involving small spherical particles, one would en- 
counter the equation 



v-[kMv*]=o. 



(12) 



By applying the method of separating variables in a 
spherical coordinate system, one finds that \f/ can be 
expressed in terms of the following characteristic 
functions, i.e., 



t=W n (r)P%(cos 8) ™ m<t>, (13) 



where W n (r) satisfies the equation 

jUjM'O ^r~\-n(n+l) K (r)W n =0, (14) 

and P^(cos 0) denotes the associated Legendre func- 
tions. The solutions to (14) for several profiles of 
k(t) are given below. 

3.1. Luneberg Lens: K(r)=2—r 2 

In this case the radial function W n (r) is expressible 
in terms of the hypergeometric function. The func- 
tion which is regular at r=Q is given by 



W n =(ff F(a,b,c; r ~) 



(15) 



where F(a,b,c;z) denotes the hypergeometric func- 
tion, and 

i 3 



3.2. Maxwell Fisheyes: /c(r) = 



(l+r 2 ) 2 



The radial function, again, is expressible in terms 
of the hypergeometric function. The function which 
is regular at r =0 is 



W n tr) = (-r^F(a,b,c;-r*), 



(16) 
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where 



n 3 . 1 






,3 

c=n +2 



3.3. Bilinear Distribution: k(t)=K. z 



r+r x 
r+r 2 



By letting ir„ = 
the equation 



-r n Z ni one finds that Z n satisfies 



Z " + L r +r+r+r+rJ Z " 

W 1 1 ) Zn =0. (17) 

In contrast to the wave ('unctions for the same pro- 
file of k(/*), it can be verified that the function Z n (r) 
has a regular singularity at infinity, hence it is of 
the type [0,4,0]. The function is, therefore, similar 
to the hyp ergeome trie function which is designated 
by [0,3,0]. It is obvious from (17) that when r 2 =0, 
corresponding to the conical lens distribution of 
Luneberg, Z n reduces to the hypergeometric func- 
tion. 

4. Equation v[^-V^]=— p/e and Its Solution 

The foundation of the electromagnetic theory deal- 
ing with a harmonically varying field in a magneto- 
ionic medium is based upon the following Maxwell 
equations: 



VXE=-j» coH, 



VXH=J+jue Q K.E, 



(18) 



(19) 



where J denotes the current density function not of 

magneto-ionic origin. The dielectric tensor k, for the 
case that the static magnetic field is pointed in the 
^-direction, is given by 



(20) 



Numerous authors have recently investigated the 
complete wave solutions resulting from (18) and 
(19). We would like to mention particularly the 
works of Bunion [1957] and Arbel [1960], where 
both use the Fourier-transform method in deriving 
the formal solution to the problem. The resultant 
integrals can be evaluated by the asymptotic method 



Kll 


*12 


' 


Kl2 


Kll 











K 33- 



to yield the far-zone field of a radiating element 
No explicit solution is yet available for the near- 
zone field except the brief description given by 
Mitt rn and Deschamps [1962] in the Abstract of the 
recent Copenhagen Electromagnetic Symposium 
(July 1962). The interesting problem of the scatter- 
ing of an electromagnetic wave by an obstacle in 
an anisotropic medium is yet to be solved. The 
problem of the scattering of a plane wave by a small 
spherical anisotropic body in air has successfully 
been solved by Berk and Lengyl [1955]. With this 
brief introduction it seems of some interest to un- 
cover whatever information we can gather from 
(18) and (19), even in a restricted form. We shall 
consider here the quasi-static solution to the electric 
field by neglecting the right term of (18), and hence 
start with 



vx#=o, 



(21) 



and the divergence equation resulting from (19) and 
the equation of continuity, i.e., 



In view of (21), one may introduce 
function \p such that 

Substituting (23) into (22) one obtains 



v.| 



6(1 



(22) 
potential 

(2.3) 
(24) 



In a cartesian coordinate 4 system, (24) can be written 
as 

Kll \w + w) +3i W = -7,>' (25) 

where we have already made use of (20). To facili- 
tate the analysis, we shall first consider the case 
where p represents that of a single charge located 
at the origin; i.e., we let 



p=qd(r-Q), 



(26) 



where 8(r— 0) denotes the three-dimensional delta 
function. Equation (25) can then be written as 

^+^+«>^- JL Hr-0). (27) 

dx 2 dy 2 dz 2 e «n 

By means of the Fourier-transform method or by 
transforming (27) into the standard form of Poisson's 
equation by letting x f = ax and y' = ay, one finds that 
the solution to (27) is given by 



(28) 



4:T6 {) K u I\ a 
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where 



E a =,[ a 2 {x 2 +y 2 ) + z 2 ]\ 



Based upon this elementary solution, the field of a 
dipole situated at the origin and of moment p pointed 
in the ^-direction can be derived. The result is 
given by 



p-R 



where 



p 4:ire K n Rl 

E=^J x 2 +y 2 +z 2 . 



(29) 



The corresponding electric field can be found by 

means of (23). The dependence of E upon the 
modified radial distance R a concurs with the principal 
term of the near-zone field of an oscillating dipole 
in a magneto-ionic medium as obtained by Mittra 
and Deschamps through the wave solution. Equa- 
tion (28) may also be used to find the characteristic 
functions pertaining to the homogeneous equation 



■ 2^ 



dx 2+ by 2 ~ 



bz 2 



(30) 



through the technique of multipole expansion. It 
is, however, more convenient to find these solutions 
by changing (30) into the standard form of the 
Laplace equation 



where 



dV d^ dv 
i)xV dyV dz 2 

X a = aX, y a = a y. 



(31) 



The characteristic functions, when expressed in a 
spherical coordinate system, are then given by 



^ nm e o HAnR n a+B n R- (n+1 ^P^(cose a )S?n^, (32) 



where 



R a =[a 2 (x 2 +y 2 ) + z 2 ] l/2 

= R[a 2 +(l-a 2 )cOS 2 d] 1/2 



cos 6 a =—-= 



cos 6 



R a [a 2 +(l-a 2 ) COS 2 ^ 1 ' 



The dipole field as given by (29) corresponds to the 
function 

In the work of Berk and Lengyl [1955], only the 
characteristic functions 



RaPKcos a )£ 4>=<*RPl(cos B)£ 



■*««{; 



are needed to satisfy the required boundary condi- 
tions for the corresponding magnetostatic problem. 
There seems to be no simple way to find the solution 
for the complementary problem of a small spherical 
dielectric body placed in an anisotropic medium. 
The only problem of this nature that can be solved 
neatly consists of a perfectly conducting spheroid 
with a contour defined by 



a 2 (x 2 +y 2 ) + z 2 =K 2 . 



(33) 



When such a spheroid is placed in an anisotropic 
medium, then an incident field defined by 



^t=— E x 

would produce a scattered field given bj 

K 3 , 
\// s =E pg Pi (cos 6a) cos 0. 



(34) 
(35) 



The author gratefully acknowledges the discussions 
which he had with Dr. H. Unz on the problem 
involving an anisotropic medium. 
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